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Abstract. We describe a practical method for finding an L— function without first finding 
the associated underlying object. The procedure involves using the Euler product and the 
approximate functional equation in a new way. No use is made of the functional equation 
of twists of the L— function. The method is used to find a large number of Maass forms on 
SL(3, Z) and to give the first examples of Maass forms of higher level on GL(3), and on 
GL[A) and Sp{4). 



1. Introduction 

We describe a general approach to finding an L— function given only a limited amount of 
information about its functional equation and Euler product. We illustrate the method by 
locating L— functions associated to Maass forms on SL(3, Z) and some of its subgroups, and 
on Sp(A,Z) and SX(4,Z). 

1.1. Summary of the method and results. We describe an approach to solving the 
following problem: find all L— functions satisfying a given functional equation. The precise 
assumptions, which also involve an Euler product and a Ramanujan bound on the coefficients, 



are given in Section |3Tj Our results show that, at least for L— functions of degree d < 4, this 
problem can be solved in a practical way, without an assumption of a functional equation 
for twists of the L— function and without first finding an arithmetic object that gives rise to 
the L- function. 

Our approach involves a new way of extracting information from the smoothed approx- 
imate functional equation for an L— function. This formula, given in Theorem 4.1, is a 
standard tool in analytic number theory. The new ingredient makes use of the fact that a 
variety of test functions can appear in the formula. This allows us to create a relatively 
small system of equations which the Dirichlet coefficients of the unknown L— function must 
satisfy. Solving that system provides the missing information about the L— functions with 
the given functional equation. 

Using this method, we have found more than 2000 L— functions of Maass forms on SL(3, Z) 
and several dozen L— functions of Maass forms on congruence subgroups of SL(3, Z). The 
method also shows the non-existence of small eigenvalues for Maass forms on SL(3, Z) in a 
manner that is more effective than previous results using the explicit formula [20] . 

We have also found more than 200 L— functions of Maass forms on SL{4, Z) and several 
dozen for Sp(4, Z). The smallest Sp(4, Z) example has the surprising property that its first 
zero on the critical line has a larger imaginary part than the first zero of the Riemann zeta- 
function [3]. We have also used this method to find, with proof, the smallest conductor of a 
hyperelliptic curve [XT] . 



The authors thank Ce Bian, Andrew Booker, Brian Conrey, Stephen D. Miller, Michael Rubinstein, Ralf 
Schmidt, and Nicolas Templier for helpful discussions. 
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Several independent tests are applied to our examples, and confirm our results in all 
cases. This gives us confidence to claim that our approximate functional equation method 
does indeed find L— functions without the need to first identify the underlying object. 

Our data will be made available at http://www.LMFDB.Org/L/degree3 and /degree4. 



2. Prior methods 

We describe methods which have been used to find automorphic forms on GL(2), as well 
as recent approaches for higher rank groups. 

2.1. Holomorphic modular forms and Fourier series methods. There are two types of 
modular forms on a Hecke congruence group Tq(N) C SL(2, Z): holomorphic modular forms 
and Maass forms. Constructing holomorphic modular forms is relatively straightforward. If 
the level, N, is small, then one can find an explicit basis of the ring of modular forms 
on Fq(N). For any level N and any weight k, one can use modular symbols [27] to produce 
a spanning set for the space Sk(T (N)). Diagonalizing with respect to the Hecke operators 
and the Atkin-Lehner operators gives the basis of newforms, which have algebraic integer 
Fourier coefficients. This process has been completely automated |24j . 

For Maass forms, the situation is quite different. Except for a thin set arising from qua- 
dratic fields, there is no known explicit construction of these functions. Thus, for most 
Maass forms one must rely on computer calculations to determine numerical approxima- 
tions [121 H3 I2H] • Recently, such computations have been proven to be correct [5] . 

The methods used to find Maass forms on GL(2) make use of the Fourier expansion and 
the transformation properties under elements of the group. For GL(d), d > 3, the Maass 
forms are functions of \{d 2 + d — 2) real variables. Their Fourier expansions [HI E] are 
(d — l)-fold sums over the integers, with an embedded sum over cosets of SL(d — 1, Z). The 
complicated nature of the sum makes it difficult to directly find Maass forms on SX(3,Z), 
although partial success has been reported by Bian and Mezhericher [TJ HH]. For SL(A, Z) 
and higher, that approach is probably infeasible. 

2.2. L— function methods. A natural alternative to dealing with the Maass form is to 
work directly with its L— function. There is no loss of information passing between those 
functions since the spectral data and Fourier coefficients of the Maass form can be recovered 
from the L— function. This approach was used successfully by Bian [HE], who made use of the 
functional equations for twists of the L— function, invoking the GL(3) converse theorem. The 
method involved solving a system of linear equations with approximately 10,000 unknowns. 
More details on Bian's method are given in Section [6j Another approach was suggested by 
Miller [12], based on the Voronoi summation formula. 

In this paper we present an approach which also involves finding the L— function, but we 
assume only one functional equation and do not make the assumption of a functional equation 
for twists. Instead, we make use of the Euler product to produce additional relations. This 
approach enables us to work with a very small system of equations, typically fewer than 30 
unknowns, although it is a non-linear system. 

In the next section we describe the results we have obtained with our method, and in 
Section [4] we give details about our algorithm. 
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3. Notation and Results 

There have been several axiomatic definitions of "L— function," most notably that of 
Selberg [25J who conjectured that the functions in the Selberg class coincide with those arising 
from automorphic representations. Below, we give a more restrictive set of axioms, which 
describes the properties that are known or conjectured to hold for L— functions associated 
to a unitary cuspidal automorphic representation of GL(d). 

3.1. L— function axioms. We consider L— functions given by a Dirichlet series 

oo 

Z / n S 

n=l 

where a n <C n e for any e > 0. The estimate on the Dirichlet coefficients is known as the 
Ramanujan bound. 



The Ramanujan bound implies that the Dirichlet series (3.1) converges absolutely for 
a = 3?(s) > 1. We assume that L(s) continues to an entire function. (In general an L- 
function can have a pole at s = 1, but the L-functions we consider will be entire.) We 
assume that the completed L— function, A(s), is entire and bounded in vertical strips, and 
satisfies a functional equation of the form 

A{s) =N%Y[ T R (s + N ) J] T C ( S + v 3 ) ■ L(s) 
i=i j=l 
= eA(l - s), (3.2) 

where iV is a positive integer called the level, and |e| = 1. Here Tk(s) = 7r _s ' 2 r(s/2) 
and r"c(s) = 2(27r)~ s r(s) where T is the Euler T-function, and A(z) = A(z). The integer 
d — d\ + 2d 2 is called the degree of the L— function. The analogue of the Selberg eigenvalue 
conjecture is that ^t(fij), 3ft (yj) > 0. 

We also assume an Euler product of the form 

L(s) = l[L p (s)l[L p (s), (3.3) 

p\N p\N 

where L p (s) = / p (p _s ) _1 with f p a polynomial with f p (0) = 1. If p \ N then we assume f p 
has degree d and satisfies the self-reciprocal condition f p (z) = x{p){~ ^) d z d f p (z~ l ), where 
X is a Dirichlet character mod N, known as the central character of the L-function. If p\N 
then we assume f p has degree < d — 1. 

An equivalent formulation of the Ramanujan bound is that if p \ N then all roots of f p 
lie on the unit circle, and if p\N then all roots of f p lie on or outside the unit circle. This 
implies a more precise form of the Ramanujan bound: 

WP"^^- 1 ). (3.4) 

3.2. L— functions of Maass forms. If F is a Maass form on SL(d,Z) with Fourier coeffi- 
cients (a n i,...,n d _i), then its L— function L(s,F) has Dirichlet coefficients a n = a^i,...,].. The 



L— function has an Euler product of the form (3.3) with u p = 1. The Ramanujan bound 



has not been proven in this case. The L— function satisfies a functional equation of the form 



(3.2) with N = 1, e = 1, d\ = d, d 2 = 0, and {fix, . . . , /id-i} = {i^i, ■ ■ ■ , i^d-i}> where the Aj 
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have simple expressions in terms of the eigenvalues of the Maass form under the generators 
of the ring of invariant differential operators, and Hd = — Ylj<d^3- The Selberg eigenvalue 
conjecture asserts that the Xj are real. This was proven by Miller [20] when N = 1 and 
d = 3. See [HI E] for more details about these L— functions. 



3.3. Results for SX(3, Z). We continue the notation of Section 3.2 with d = 3. 

By permuting the Xj and possibly replacing the L— function by its dual (conjugate), we 
can assume Ai > A 2 > and A 3 = — Ai — A 2 . Thus, we specify the functional equation of the 
L— function, equivalently the eigenvalues of the Maass form, by a point (Ai,A2) below the 
diagonal of the first quadrant of IR 2 . The symmetric squares of SL(2, Z) Maass forms have 
A 2 = 0. 

Weyl's law for SL(3, Z), proved by Miller [20J, asserts that 

#{(A 1; A 2 ) : Xj + AiA 2 + X\ < T} ~ v6L W y | (3 5) 

(47r)3r(|) 

where X = SL(3, Z)\SL{3, M)/SO(3, M). Miller also proved the Selberg eigenvalue conjec- 
ture for SL(3, Z) in the form A 2 + AiA 2 + A 2 . > 80. (Miller proved slightly more, which is 



described in Figure 3.1 below 



The purpose of this paper is to describe a new method for finding L— functions satisfying 



the axioms given in Section 3A_ Our results fall into two categories: an extension of Miller's 
result that certain regions do not contain the eigenvalues of any SL(3, Z) Maass forms, and 
the experimental determination of a large number of L— functions that appear to arise from 



SL(3, Z) Maass forms. Both types of results are displayed in Figure 3.1 



The dark-shaded region in Figure 3.1| was shown by Miller [20], using the explicit formula, 
to contain no eigenvalues (Ai,A 2 ) of any SL(3, Z) Maass forms. We use the approximate 
functional equation to prove a stronger result. 

Theorem 3.1. Assume the Ramanujan bound \a p \ < 3. There are no L— functions with 
functional equation ( 3.2[ ) with d = 3 and N = 1, therefore no SL(3, Z) Maass forms, with 



spectral parameters (Ai, A 2 ) in the lightly shaded region shown in Figure 3.1 



The proof, which is a new application of the approximate functional equation, is given in 
Section [7l 



Also shown in Figure 3.1 are the 124 lowest (in lexicographical order) pairs (Ai,A 2 ) of 
L— functions with functional equation and Euler product as described in Section |3.2| with 
d = 3, which were found in our search. We believe that the figure shows all cases with 
Ai < 20, and may include all with Ai < 24. Numerical values for the first 15 examples, and 



their first two Dirichlet coefficients, are listed in Table 3.1 

Ce Bian (personal communication) has confirmed several dozen of our examples and com- 
puted a few thousand Dirichlet coefficients in each of those cases. 

We have found more than 2000 spectral parameters (Ai,A 2 ) and associated Dirichlet co- 
efficients. These will be made available at |http://www.LMFDB.org /L/deg ree3| 

3.4. Subgroups of SL(3, Z). We applied the approximate functional equation method to 
find L— functions associated to Maass forms on subgroups of SX(3,Z). In terms of the 
functional equation, we have d\ = 3, <i 2 = 0, and N > 1. For p\N, the local factor in the 



Euler product (3.3) will have degree at most 2. In our calculations, we assume an arbitrary 
degree-2 local factor at the primes dividing the level: L p (s) = (1 + A p p~ s + B p p~ 2s )~ 1 . It 
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Figure 3.1. The axes are Ai and A 2 and the dots indicate (Ai,A 2 ) for the 
parameters in the functional equation of SL(3,Z) L— functions, for the first few 
examples we found. The darkest region is Miller's excluded region |20J, and we 
prove, assuming the Ramanujan conjecture, that the lighter shaded region contains 
no eigenvalues. 



Ai 


A 2 


a 2 


a 3 


13.59658451 
14.14163558 
15.31863407 


4.76468206 
2.38038848 
4.19173391 


1.04846245 - 0.37523963?' 
-0.10524097 + 0.75072694?' 
0.33541208 + 0.34590083? 


-0.49904094 - 0.27897508?' 

1.23599391 - 0.03911217? 
-0.27260607 - 0.97283444?' 


15.74069912 
16.05436164 
16.40312474 


7.85232504 
1.98365457 
0.17112189 


-0.39016397 - 0.62517821?' 
-0.85221921 + 0.418004102 
-0.42168648 + 1.06796797? 


0.75666087 - 0.16479552?' 
0.67428665 - 0.33153749?' 
-0.76802216 - 1.31329241?' 


16.75957309 
17.34401833 
17.42523780 


4.03941275 
7.13419625 
5.86543671 


0.43484942 - 0.64408778? 
0.55207716 - 0.70388045? 
-0.86407728 - 1.64923555? 


0.67804925 - 0.27020173? 
0.03833728 + 0.25697818? 
-0.22106187- 0.90868355? 


17.50092882 
17.65047391 
17.86523805 


1.90552951 
7.30146768 
10.34945597 


-0.28445287 + 0.01165958? 
-0.70569544 + 0.37307990? 
-0.86106525 + 0.35416339? 


0.42424909 + 0.46219014? 
-0.29251139 + 1.07644911?' 
-0.13894817 + 0.24569291? 


18.06174615 
18.08927092 
18.19481530 


8.34053395 
5.30177649 
0.49254810 


-1.15143879 - 1.75641008? 

0.32185117 + 0.32366378? 
-1.15442278 + 1.90981460?' 


-0.49129679 + 0.21151351? 
-0.36219655 - 0.00596403? 
0.33665482 - 0.10748733?' 



Table 3.1. Spectral parameters (Ai, A 2 ) and sample Dirichlet series coefficients 
for the first 15 L— functions of SX(3, Z) Maass forms. The numbers are truncated 
to 8 correct decimals. 
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may happen that B p and/or A p is zero, if the local factor had degree or 1. By making no 
assumptions on the local factors, we have the opportunity for an independent check on our 
calculations. 

We applied the approximate functional equation method to search for L— functions of 
levels N = 2 and 3, and trivial central character, and our search found nothing. This seemed 
at first to be problematic because when doing numerical experiments it can be difficult to 
distinguish a true non-result from a non-result due to a error in the computer program. In 
this particular case, it turns out that we experimentally discovered the theorem, which is 
relatively straightforward but does not seem to be in the literature: there are no Maass forms 
on GL(3) which give rise to L— functions of squarefree level and trivial central character. 
The most natural small index subgroups of SL(3, Z) have outer automorphisms (analogous 
to the Fricke involution for To(iV)) which give level p 2 in the functional equation. 

So we searched for L— functions of levels N = 4, and 9 and found examples in each 
case. For square level, the corresponding subgroups have outer automorphisms of order 3. 
Therefore, the signs which appear in the functional equation should be 3rd roots of 1, and 
we found examples of all three signs. For level p 2 , in every case the local factor we found 
at p was L p (s) = (1 — ep~ l ~ s )~ l , where e is the sign of the functional equation. Note that 
we assumed a completely general degree 2 local factor, but then found that it actually had 
degree 1. 

Ralf Schmidt (personal communication) informed us that this is the local factor from the 
Steinberg representation for GL(3), which has level p 2 . This is a strong independent check 
of the validity of our calculations. 

The first few L— functions for level 4 are given in Table 3.2 The expected lifts from GL(2) 
were found in our search. 



3.5. Results for Sp{4, Z) and GL(4, Z). We have also used our method to find degree 4 
L— functions, that is, L— functions with functional equation (3.2) and Euler product (3.3) 
with d = 4. At present, we consider the case of level N = 1; we will address the situation of 
higher level in a subsequent paper. 

It is natural to distinguish the special case where Ai = —A3 and A2 = — A4, and the 
Dirichlet coefficients a n are real. We refer to these as u Sp(4, Z) L— functions" because we 
believe they are the spin L— functions of the analogue of Maass forms on Sp(4, Z). Since the 
functional equation depends on two real parameters, we can provide a plot of our results, 
analogous to the case of SL(3, Z). Figure 3.2 shows the first 58 examples that we have found 
for Sp(4,Z) Maass form L— functions. The first 15 of these are also listed in Table 3.3 with 
spectral parameters and the first four prime coefficients. 

In Figure 3^ one observes that the some of the first few L— functions lie approximately 
along certain lines. A similar phenomenon can be seen in Figure 3.1[ 

The first Sp(A, Z) L— function, which has (Ai,A2) = (12.46875,4.72095) has the curious 
property that its lowest zero on the critical line is at | ± 14.496i, which is higher than the 
lowest zero of the Riemann zeta-function. This is discussed further by Bober, et al. j3]. As a 
check on these calculations, we also compute the standard (degree 5) L— function using the 
data from the spin (degree 4) L— function (see section [6]). 
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e 


Ai 


A 2 


a 2 


a 3 


1 


8.23979752 


2.64122672 


1/2 + 5 


1.06604225 + 0.48105885i 


1 


10.64614640 


5.40948610 


1/2 + 5 


0.06047305 - 0.43014710i 


1 


11.72327411 


2.02460527 


1/2 + 5 


0.59406181 - 0.2755031H 


1 


12.42065536 


4.72210152 


1/2 + 5 


0.99900525 - 0.54164255i 


/ 1 \ 

e (§) 


9.63244452 


1.37406028 


e(l)/2 + 5 


0.15012282 + 0.78657327i 


e (§) 


10.69551431 


3.18301083 


e(l)/2 + 5 


0.62468273 - 0.09470423i 


/ 1 \ 

e (|) 


11.25849528 


1.13707825 


e(l)/2 + 5 


-0.11242070 + 0.10083473i 


/ 1 \ 


12.64900811 


0.92183852 


e(l)/2 + 5 


-0.80493907 - 0.12717437i 


<l) 


8.95466251 


2.93659153 


e{\)/2 + 5 


0.36025530 - 0.35631998i 


<D 


10.06466234 


4.51670022 


e{\)/2 + 5 


-0.90149321 - 0.57545547i 


<D 


10.73782624 


2.37415775 


e(l)/2 + 5 


-0.53488262 + 0.88247158i 


<l) 


10.82653286 


0.77122693 


e(l)/2 + 5 


0.58740526 - 0.33415990i 



Table 3.2. Spectral parameters (Ai,A 2 ), sign (e), and sample Dirichlet series 
coefficients for L— functions of level 4 of degree 3 Maass forms. Here e(x) = 
exp(27rix) and \5\ < 10 -9 . The numbers are truncated to what we believe are 
8 correct decimals. 
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20 



Figure 3.2. The axes are Ai and \ 2 and the dots indicate (Ai, A 2 , — Ai, — A 2 ) 
for the first few Sp(A,Z) L— functions. Note that we may assume Ai > A 2 . 

We have also found more than 200 L— functions of SX(4, Z) Maass forms that are not 
Sp(A, Z) L— functions. We list the eigenvalues and first few coefficients of the smallest ex- 
ample. We believe that all the given digits are correct: 

(Ai,A 2 ,A 3 ) = (13.048202385582,6.726210585782,-4.367763225255) 

a 2 = - 0.1912480688849 - 0.587500805369 i 

a 3 = - 0.652904479858 + 0.96315195085 i 

a A = - 0.5309393169 + 0.22471678899 i 

a 5 = - 0.3125111187 -0.0112785467i. 
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Ai 


A 2 [ 


a 2 


03 


a 5 


a 7 


12.46875226 
14.09372017 
14.72547970 


4.72095103 
8.83036460 
4.02669357 


1.34260324 
-0.31644276 
0.73680321 


-0.18745190 
0.05294768 
0.99830784 


-0.0016279 
0.0261060 
0.2677819 


0.228229 
0.952439 
-0.031085 


14.83756247 
15.23336181 
15.62467181 


6.41243964 
8.25653712 
9.46217647 


0.18715348 
0.88602292 
-1.72583459 


-0.21546510 
1.76959583 
1.11618173 


0.5492948 
-0.1217794 
-0.6444860 


0.582942 
-0.179071 
-0.009432 


15.63453910 
15.96343790 
16.00467179 


5.57458793 
11.59606479 
8.12341090 


2.30065434 
0.17708238 
-0.59748302 


0.16113045 
-0.36379926 
-1.00280803 


-0.4928333 
-1.3728032 
-1.5990447 


-0.410839 
0.233579 
0.398796 


16.17300694 
16.19395973 
16.54574513 


3.71016260 
6.05889766 
9.57822088 


0.60584003 
0.04742685 
0.27348521 


0.74786457 
0.01653540 
0.97263132 


-1.1437477 
1.6223066 
1.4040048 


0.469017 
-0.080626 
-0.955692 


16.77928785 
16.78915308 
16.81031209 


10.78021484 
6.16376335 
9.58414788 


-1.50813572 
-1.29328034 
0.22614271 


-0.48797568 
-1.51751705 
-0.89126171 


0.1694122 
-0.0818331 
0.0562885 


-0.471615 
-0.784244 
-0.906897 



Table 3.3. Spectral parameters (Ai, A 2 ) and sample Dirichlet series coefficients 
for 15 L— functions of Sp(A, Z) Maass forms. The numbers are truncated to what 
we believe are 8 correct decimals. 



Numerical values of the spectral parameters and a 2 for the first 15 examples we found 



are listed in Table 3.4 They are normalized so that Ai > A 2 > A3 > A4 and |Ai| < | A4 1 
and in lexicographical order. We don't claim that this list necessarily contains the first 15 
L— functions with this ordering. The full list, including the Dirichlet coefficients, will be 
made available at |http://www.LMFDB.org /L/de gree4| 

3.6. Other types of degree 4 L— functions. As a test of the approximate functional 
equation method, we now apply it to a different type of degree 4 L— function. In this 
example, we search for L— functions satisfying a functional equation of the form 

AO) = T c (s + l)T c {s + k— |)L(s) = A(l - s). (3.6) 

That is the functional equation of the spin L— function of a Siegel modular form of weight k on 
Sp(4, Z). As a demonstration of the method, let k = 20. It is known that there is exactly one 
such L— function, associated to the cusp form denoted T 2 o- Its first few Dirichlet coefficients, 
in the arithmetic normalization, are the following rational integers [26J : 

2 37 / 2 a 2 = - 840960 
3 37 /2 as = 346935960 

4 37 / 2 a 4 = 316975677440 

5 37 / 2 a 5 = _ 5232247240500. (3.7) 
We applied the method of Section|4|to search for the Dirichlet coefficients of an L— function 



satisfying the functional equation (3.6) with k = 20, and having a degree-4 Euler product. 
We used the first 90 terms of the Dirichlet series, which is 30 unknowns because a degree-4 
local factor is determined by a p and a p 2. We chose 30 triples (| + it, ^1,^2) with < t < 33, 
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Ai 


A 2 


A 3 


a 2 


13.04820238 
13.74207582 
14.94493797 


6.72621058 
8.91090555 
10.17319896 


-4.36776322 
-6.67116480 
-7.63795646 


-0.19124806 - 0.58750080i 
-0.38362321 + 0.8827767H 
0.50916811 + 0.08060143i 


15.82354814 
15.84333613 
15.94944916 


5.98442035 
8.14020386 
11.22357307 


-3.32054209 
-6.39808715 
-7.46849155 


0.06707773 - 2.25215645i 
-0.93727868 + 0.45045257i 
-0.01250803 - 0.16785927i 


15.99896071 
16.04868870 
16.17564264 


5.55250290 
7.44354895 
9.31222108 


-4.28979199 
-5.74749580 
-7.99494686 


0.76641492 + 0.72684214i 
0.39710090 - 1.01290587? 
-0.25094605 - 1.40588826? 


16.72584930 
16.89972715 
16.92627200 


10.57549478 
2.27258771 
12.47850050 


-8.05016487 
-6.03583588 
-10.63475738 


-1.81300469 - 0.42294306i 
0.55659019 + 0.92818484i 
0.33041049 + 0.59891916? 


16.95696382 
17.16046064 
17.27511800 


9.53674609 
5.51936319 
7.18638875 


-7.32770475 
-3.67883841 
-4.38401627 


-0.66772182 + 0.75485571? 
-0.31914164 - 0.69706013? 
-0.13829583 + 1.09540674?' 



Table 3.4. Spectral parameters (Ai, A 2 , A 3 ) and a 2 for L— functions of SL{A, Z) 
Maass forms. The numbers are truncated to what we believe are 8 correct decimals. 



selected so that the truncation error in eliminating a n for n > 90 is less than 10 -12 . After 
using the Euler relations to convert to a non-linear system, we solved the system 1000 times 
using the secant method with random starting values. Only one solution was found (and it 
was found repeatedly), which had: 

2 37/2 a 2 = - 840960.000014382993 

3 3? / 2 a 3 = 346935960.202402521 

4 3? / 2 a 4 = 316975677305.428055 

5 37 / 2 a 5 = - 5232247227434.34870. (3.8) 

We see that the (rescaled) a 2 is convincingly close to an integer, and it is the correct integer. 
Replacing a 2 by its apparent exact value, and the re-solving for the other coefficients, we 
find: 

3 37 / 2 a 3 = 346935959.999827467 
4 3? / 2 a 4 = 316975677440.106735 

5 37/2 a 5 = - 5232247240497.56648. (3.9) 

Now the rescaled a 3 is convincingly close to an integer, and again it is the correct integer. 
Replacing a 3 by its apparent exact value, and solving again we find: 

4 3? / 2 a 4 = 316975677439.993180 

5 3? / 2 a 5 = - 5232247240498.27543. (3.10) 

This time, the situation is less clear. The rescaled a 4 certainly appears to be close to an 
integer; in fact the correct integer. But we had truncated the Dirichlet series to give 12 digits 
of accuracy. Since the true value of 4 37//2 a 4 is an integer with 12 digits, our calculation should 
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not have been adequate to determine that integer exactly. It is somewhat surprising that we 
obtain a few more correct digits than we should have expected. 

This example illustrates that the approximate functional equation method should be ap- 
plicable to other types of L— functions. We will address these extensions of the method in a 
subsequent paper. 



4. The method: approximate functional equations 

Our method of locating L— functions is based on a smoothed approximate functional equa- 
tion. The approximate functional equation is a tool that typically is used to prove results 
about L— functions (such as moments), or to calculate L— functions [22]. Here we turn it into 
a method of locating L— functions by exploiting the fact that the approximate functional 
equation contains a free parameter. This approach was inspired by Michael Rubinstein's 
Lcalc program [22, 23J. 



4.1. Smoothed approximate functional equation. Let 

oo 

^— ' n s 

n=l 

be a Dirichlet series that converges absolutely in a half plane, dt(s) > <J\. 
Let 

A(s) = Q s (f[r(K j s + \ j )jL(s), (4.2) 

with Q, Kj G M + , 3?(Aj) > 0. Note that the functional equation we present here is somewhat 
more general that given in the previous section. We assume that: 

(1) A(s) has a meromorphic continuation to all of C with simple poles at S\, . . . , St and 
corresponding residues r±, . . . , Tg. 

(2) A(s) = eA(l - s) for some e G C, |e| = 1. 

(3) For any <r 2 < 03, L(a + it) = 0(expt A ) for some A > 0, as \t\ — > 00, a 2 < a < cr 3 , 
with A and the constant in the 'Oh' notation depending on 02 and 03. 

To obtain a smoothed approximate functional equation with desirable properties, Rubin- 
stein [22] introduces an auxiliary function. Let g : C — > C be an entire function that, for 
fixed s, satisfies 

\A(z + s)g(z + s)z~ l \ -)■ 

as 1^(^)1 — > 00, in vertical strips, — xq < < Xq. The smoothed approximate functional 
equation has the following form. 

Theorem 4.1. For s ^ {s±, . . . , si}, and L(s), g(s) as above, 

l I N OO OO 

L — ' s — Sk L — ' n s L — ' n 1 s 

k=l n=l n=l 
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where 



fi(s,n) 



/ 2 (1 - s,n) 



1 
1 



/ II r M-- + s) + A i )z~ 1 ^(s + z)(Q/n) z dz 

a 

U r(^(^ + 1 - s) + Aj> _1 #(s - ^)(g/n) 2 dz 



(4.4) 



u>zt/i z/ > max {0, — 3?(Ai/ki + s) 



-3?(A a /K a + s)}. 



Note that, despite its name, the approximate functional equation is an exact formula for 
the L— function throughout the complex plane. 

We assume L(s) continues to an entire function, so the first sum in (4.3) does not appear. 
For fixed Q, k, A, e, sequence {a n }, and function g(s), the right side of (4.3) can be evaluated 
to high precision. In the next section we describe how we use this observation to model the 
functional equation. 

4.2. Modeling the functional equation. The Z-function (or Hardy Z-function) of an 
L— function is an analytic function which is smooth and real on the |-line and satisfies 
\Z{\ + it)\ = \L{\ + it)\ for t £ R. In terms of the completed L— function (4.2), 

a 

Z{\ + it) = (eQ)-i Y[ irfod + it) + Aj)!" 1 ■ A(i + it), (4.5) 

3=1 

when tel. That formula is not a valid definition of Z(| +ii) if t ^ R, but it is sufficient for 
our purposes because we will only evaluate the L— function on the critical line. We will find 
the Z-function to be convenient because all calculations will involve real numbers instead of 
complex numbers. 

Fix all the parameters in the functional equation of an L— function and write (4.3) as 
A(s)g(s) = AppFE(s, Q, k, A, e, {a n }, g), where we suppress the contribution of the poles 
because we will assume our L— functions are entire. Use (4.5) to rewrite (4.3) as 

a 

Z(s) = Z(s, g) = g{s)- l {eQfz J] \T(K jS + ) | 1 A PP FE(s, Q, «, A, e, {a n }, g) 



g(s) 1 RHS(s,Q,k, A,e, {a n },g), 



(4.6) 



say. 

The only free parameters in the above equation are the complex number s and the func- 
tion g. The left side of (4.6) does not depend on the function g. Thus, for each fixed s, if 
we evaluate the right side of (4.6) with two different g functions gi, g2, then setting them 
equal gives an equation in the Dirichlet coefficients of the L— function. 

Thus, if we fix So = | + Hq and choose two functions gi, g<i then 

gi^r 1 RHS(s , Q, k, A, e, {a n }, g x ) - g 2 (s o y 1 RHS(s , Q, «, A, e, {a n }, g 2 ) = (4.7) 

is a linear equation which the Dirichlet coefficients of the L— function must satisfy. 

We can produce many equations by choosing various (s, pi,^) combinations. If we have 
chosen functional equation parameters which correspond to a genuine L— function, then its 
Dirichlet coefficients will be a solution to that system. If we have chosen parameters for which 
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there is no L— function, then it may be possible to detect this through an inconsistency in 
the system. 

We expand on this approach in the next section. 

4.3. The method. It is easiest to first think of the functional equation as given, and focus 
on determining the Dirichlet coefficients. For the case of Maass forms, we do not actually 
know the functional equation in advance: we must find the spectral parameters Ai, . . . , A^-i 
which determine the functional equation. This is an added complication which we address 
in Section [5TT1 

Assume we are given the functional equation of an L— function and we want to find its 
Dirichlet coefficients. The main idea is to first choose various Sj = | + itj, and pairs of test 
functions and gj j2 to form the system of linear equations {Z(sj,gj t i) — Z(sj,gj t2 ) = 0}. 

Step 1. Forming the linear system 

Typically we choose g(s) = e~ ll3s+as with /3, a G R. This is an allowable test function in 
the approximate functional equation provided a > 0, or a = and \(3\ < 

We illustrate the ideas with a numerical example. The example concerns the L— function of 
a Maass form on SL(3, Z), so in Q4.2[ ) we set Q = 7r~ 3 / 2 , a = 3, K\ = k 2 = k 3 = |, and e = 1. 
This specifies all the parameters in the functional equation except Xj. We illustrate below 
with the test case A = iX 2 , 1X3) = (8.4i, 14.2i, — 22. 6i). Those parameters were chosen 
randomly, so there is no reason to expect that this corresponds to an actual L-function. In 
fact, in Section [7] we will prove that there is no L-function with that functional equation. 

Note that we must treat the real and imaginary parts of the Dirichlet coefficients as 
separate real variables, so we write a n = b\{n) + 262 (n) with b\ (n), b 2 (n) G R. 

We apply the procedure described in the previous section. This involves choosing a com- 
plex number s and two test functions g±, g 2 . Choose s = | + i, and let gi(s) = e~ 4S//4 . A 
numerical calculation finds, 

gi(so)' 1 RHS(s , Q, k, A, e, {a n }, g x ) = (4.8) 

1.2860 + 7.79076i(2) + 1.9245 63(2) + 0.6653 61 (3) - 1.0163 6 2 (3) 

+ 0.03319 h (4) - 0.1197 6z(4) + ■ ■ ■ + 4.17 • 10~ 6 61(8) - 2.37 • 10" 7 6 2 (8) + • • ■ . 

The terms decrease rapidly: the coefficients of &i(15) and 62(15) are less than 2 ■ 10 -13 . The 
exponential decay in the contributions of the Dirichlet coefficients is a general phenomenon jH 

EEB]. 

With the same Sq, now choose g 2 (s) = e~ ts l 2 . This gives 
g 2 {s )~ l RHS(s , Q, k, A, e, {a n }, g 2 ) = (4.9) 
2.02524 + 2.21176i (2) + 1.8524 6 2 (2) + 0.7408 61 (3) + 0.2210 62(3) 
+ 0.1019 6i(4) + 0.0922 6 2 (4) + • • • 2.66 • HT 6 61(8) - 0.00 00 143 6 2 (8) + • • • . 



Subtracting (4.8 ) from (4.9 ), we obtain the following equation which the Dirichlet coefficients 
of the L— function must satisfy, if such an L— function exists: 

= 0.7392 - 5.5790 6 X (2) - 0.07216 6 2 (2) + 0.07541 61 (3) + 1.23745 6 2 (3) (4.10) 

0.06876 6i(4) + 0.2119 6 2 (4) + 1.51 ■ HTH^S) - 0.00 00 1 41 6 2 (8) + • • • . 

We have just found one equation which the Dirichlet coefficients of the L-function (if it 
exists) must satisfy. The next steps involve truncating the system so that we are dealing with 
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finitely many unknowns, and choosing various (so ) fi , i ! fi , 2) combinations to form a system of 
equations. 

Step 2. Truncating the system 

The linear equation from Step 1 has infinitely many unknowns. However, the coefficients of 
those unknowns decrease rapidly. Thus, once we select a working precision, we can truncate 
the system with an error less than our working precision. 

For the remaining steps, we assume that a working precision, 5, has been selected (for 
example, 5 = 1CT 12 ), and a given number of Dirichlet coefficients, J, have been selected (for 
example, J = 30, so we need only consider a 1? . . . , a 30 ). 

Step 3. Choosing the parameters 

Having selected a working precision S and a number of Dirichlet coefficients J, there are 
limits to which (s,gi,g 2 ) combinations can be used. Namely, we require that the error in 
the approximate functional equation be less than 5 when we use these particular values of s 
and gj and omit those terms a n with n > J. Note that this requires the Ramanujan bound 
or some other estimate on the Dirichlet coefficients. 

In particular, given 5 and J, for each function g, there is a limited range of s (which may 
be empty) which meets the requirement of the chosen truncation error. 

In practice, often it is effective to fix gi and g 2 , and to vary Sj = \ +itj to obtain different 
linear equations. However, in order to keep the system of equations as well conditioned as 
possible, one should aim to have the same rate of convergence in all equations. This can be 
achieved by varying the functions g\ and g 2 with the Sj. 

Step 4. Using the Euler product 

Except in some low-degree cases, the linear system will have infinitely many solutions. For 
example, if there are two L— functions with the same functional equation, then clearly the 
linear system will have infinitely many solutions. 

In order to get around this, we make use of the Euler product and convert to a non-linear 
system. The difficulty of working with a non-linear system is partially offset by the fact that 
the number of unknowns will be much smaller. The Euler product immediately gives that 
all Dirichlet coefficients are determined by the coefficients at prime powers. Furthermore, 
for a degree-<i Euler product, the local factors are determined by [d/2] parameters. For 
example, for degrees 2 and 3, all coefficients are determined by a p . And for degrees 4 and 5, 
all coefficients are determined by a p and a p 2. 

We apply the multiplicative relations a 6 = 0203, au = 0304, etc., to reduce the number 
of unknowns. Then we apply the recursions arising from the shape of the local factors 
(for example, a p 2 = a? p — Op~ in the case of SX(3, Z) Maass forms) to obtain a system of 
equations with a minimal number of unknowns. So, for example, for SL(3, Z) Maass forms 
the unknowns will be {a p : p prime, p < J}. 

Note: we refer to the a n as the "unknowns," but in practice we must treat the real and 
imaginary parts of a n , denoted b\ (n) and b2{n), as separate real unknowns. Since a± — 1, 
the original linear system has 2( J — 1) real unknowns. 

Step 5. Solving the non-linear system 

The non-linear system will have many solutions. Most of these are of no interest and are 
simply accidental solutions. 

We use the secant method to solve the system. This is done repeatedly, typically between 
100 and 1000 times, with different starting points. The starting points are always chosen 
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to have values within the Ramanujan bound since this is supposedly true for the interesting 
solutions. Typically 0-5 different solutions are found. In most cases, the size of the a n in the 
solutions will grow rapidly with n. For the accidental solutions, the growth is (almost) as fast 
as the rate of decrease in the coefficients. For the solutions that are close to approximations 
of true L— functions, this is not true. Here the growth in the size of the a n is smaller. Thus, 
the Ramanujan bound is an important ingredient to the method because it allows us to 
eliminate unwanted solutions. However, there are more efficient ways to avoid extraneous 



solutions, see Section 5.1 



In practice, it turns out that it is better to begin with not too many unknowns in order 
to find any of the solutions at all. Once a solution has been found, the number of unknowns 
may be increased and a new larger system is solved with the previous solution as a starting 
point. Typically we start with 10 < J < 40 and end with 60 < J < 150. 

The above procedure describes a method for finding the Dirichlet coefficients of an L— function, 
given its functional equation and the shape of its Euler product. In the case of Maass form 
L— functions, there is the added difficulty of finding the spectral parameters which determine 
the functional equation. This is discussed in the next section. 

5. Searching for the spectral parameters 

When searching for L— functions of Maass forms on SL(d, Z) or one of its subgroups, 
the spectral parameters in the functional equation are not known in advance. The greatest 
challenge is to find good approximations to these parameters. 

5.1. Initial scan. The first step is to fix a box in which to initially search for spectral 
parameters Xj corresponding to an L— function. Typically, we use boxes of side length 2, so 
for SL(3, Z) and Sp(4, Z) this will be a square (with side 2), and for SX(4, Z), a cube. In 
this box we make a grid of equally spaced points with a typical step size of 1/5 to 1/25. 

Next, we choose the number of Dirichlet coefficients, J, to use and set the working precision 
(estimated maximal size of truncation error) to 5 = 0.5 • 10 -4 . We fix test functions g\ and 
g 2 of the form <fc(s) = e~ ll3iS+as . After these choices are made, we determine the range 
t G [Ti,T 2 ] of s = | + it on the critical line for which we have a truncation error less 
than 5 in the approximate functional equation. If this range is too small, we increase J 
until we achieve a large enough range. Typically the range [T 1 ,T 2 ] will be approximately 
[2 min{ A-,- } , 2 max{ A-,- }] . 

Now we determine the number of unknowns, m, of the non-linear systems of equations. 
This will depend on whether we assume the L— functions are self-dual (meaning that the 
Dirichlet coefficients are real), and on the degree of the local factors of the Euler product. In 
the case of SL(3, Z), m is twice the number of primes less than or equal to J. For Sp(4, Z), 
it is the number of primes or squares of primes less than or equal to J. We choose m + k 
(typically with k = 8) values of Si = | + 1 , in the range T\ < tj <T 2 . Of these points, m of 
them will be used to make the systems of m non-linear equations in m unknowns, and the 
other k will be used to detect candidates for first approximations to the spectral parameters, 
as described in the next subsection. 

Now for each point in our grid of spectral parameters Xj, we go through the following 
procedure. Since the parameters Xj and test functions g\ and g 2 are fixed, each of the points 



Si gives a linear equation of the form (4.7). The Euler product is used to turn these into 



m + k non-linear equations with m unknowns. We choose m of these equations (same for 
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all points in the grid) and then use the secant method to solve this non-linear system of 
equations. The system of equations is solved a fixed number of times (typically 50) with 
starting points chosen randomly each time. In addition to these random starting points, we 
also use solutions that have already been found at adjacent points in the grid as additional 
starting values. Sometimes the system with a given starting point will converge to a solution 
with given accuracy within a set number of steps and sometimes it will not. For some points 
in the grid it will converge every time and for some never. These solutions are considered 
"the same" if the first few coefficients agree to a certain accuracy. So finally we will have a 
number (typically 0-5) of different solutions at each points in the grid. 

5.2. Refining the search. Once we have a number of different solutions for each point in 
the grid, we can find candidates for the functional equation parameters. This is when the 
k extra equations come into play. Suppose the spectral parameters used in the system of 
equations are close to those of an L— function and the Dirichlet coefficients of the solution 
are also close to those of the L— function. Then those coefficients should also be close to 
solutions to the k extra equations. Hence we want to find spectral parameters for which 
there is a solution that makes the k extra equations close to consistent. We do this by 
substituting the solution found at each point into the extra equations, forming "indicators" 
which should measure how far the point is from a true L— function. We then interpolate to 
determine where the indicators approximately vanish. 

5.2.1. Collecting nearby solutions. Since there may be several solutions at each point, some 
of which may be extraneous, we first collect those solutions which are likely to correspond 
to the same L— function. 

Fix a point P in the grid and choose D adjacent points in the grid where D is the dimension 
of the grid. In the case of SX(3,Z), we have D = 2 and choose the points just above and 
just to the right of P. For each of the solutions at P we identify the solutions at the adjacent 
points that are closest. Here 'closest' in the sense that the first few coefficients differ by as 
little as possible. Now we have sets of related solutions at D + 1 adjacent points, and each 
set could potentially correspond to an L— function. 

5.2.2. Comparing nearby solutions. Now we use each set of solutions and the k extra equa- 
tions to form indicators, and then interpolate to determine where those indicators (approx- 
imately) vanish. 

Choose D of the k extra equations at each of the D + 1 points. Substitute the related 
solutions of a set into the D extra equations to get D complex numbers (indicators) at each 
of the D + 1 points. Using linear interpolation we find a point where the indicators are 
approximately 0, which gives us a candidate for an L— function. This is done for each of the 
(^) choices of extra equations, so we get candidates. If a large enough proportion (50%) 
of these candidates is close enough to P this is considered a good candidate. We evaluate the 
average (of the central 50%) of these candidates and consider this to be an initial candidate 
for the spectral parameters of an L— function, which is passed on the the next step of the 
process. We also save an interpolation of the first few Dirichlet coefficients of the solutions in 
the set. These coefficients will serve as a starting point when trying to improve the precision 
of this candidate. 

This process is repeated for each of the related sets of solutions at P and then for each 
of the points in the grid. The result of this process is a list of initial candidates for spectral 
parameters of L— functions with a first approximation of its Dirichlet coefficients. 
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5.3. Zooming in on a candidate. The final step is to try, for each of the initial candidates, 
to zoom in to a good approximation for the spectral parameters. This is done in steps where 
the number of coefficients is gradually increased and hence also the size of the truncation 
error decreased. Each of these steps is similar to the linear interpolation that gave the initial 
candidate. A system is created and solved at D+l adjacent points with the approximation of 
the Dirichlet coefficients as a starting point. If there are solutions at each of the points, then 
we get (^j) new approximations in the same way as before. If there is no solution or if the 
average of the new approximations is far from the first then we stop pursuing that possible 
L— function, otherwise we continue. This is repeated with increasingly smaller steps and 
more coefficients. This process is deemed a success when the maximum difference among the 
(d) approximations is less than a given bound (we typically require 8 correct decimals). The 
average (of the central 50%) of these approximations is saved as a highly probable candidate 
for spectral parameters of an L— function. 

In practice, this process converges very quickly, typically yielding one or two additional 
decimal digits in each step. We have implemented the above procedure in Mathematica and 
will make our data available online [T7] . 

6. Checks on the results 

In this section, we discuss some independent checks which support the claim that our 
procedure finds approximations to L— functions of Maass forms. 

6.1. Functional equation for twists. Bian [2] computed the first four examples of L— functions 
of Maass forms for SX(3,Z). His method uses a similar starting point, but he also uses the 
functional equation of many twists of the L— function. This has the drawback of requiring 

a large number of coefficients. Already for the first few examples, his method requires a 
linear system of equations with up to 10000 unknowns. Our search found his four examples 
and a comparison of spectral parameters and Dirichlet coefficients finds agreement to within 
the claimed precision. Our method gives a higher precision in the spectral parameters and 
the first few coefficients, but Bian's method produces more coefficients. Finally, Bian has 
verified the first 50 SL(3, Z) L— functions found in our search. 

6.2. Properties of the coefficients. Another check is that in our method we leave 032 
(and also 064) as unknowns. Hence we obtain a computed value of 032 that we can compare 
with the the value derived from the computed value of 02- For all our examples for SL(3, Z) 
these values agree to at least three decimal places and in some cases (with larger spectral 
parameter requiring more coefficients) to as many as eight decimal places. This also gives 
an independent estimate for the accuracy of the coefficients, which compares well with our 
estimate based on differences in the approximations in the computation. 

A strong check involving the coefficients is that the local factors at the primes dividing 
the level for GL(3) L— functions agree to high precision with the predicted value. This is 
discussed in Section 13.41 

6.3. Lifts. A check which we consider to be particularly strong is the analytic continuation 
and functional equation for higher degree lifts of our L— functions. 

For various GL{3) examples, we computed the symmetric square of the L— functions, 
which has degree 6. That is, we used the functional equation and Dirichlet coefficients of 
the degree 3 L— function to determine the functional equation and Dirichlet coefficients of 
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the degree 6 L— function. We then used the method in Theorem 4.2 of [13] to verify that the 
L— function satisfies the conjectured functional equation. In all cases we find that degree 6 
L— function satisfies the expected functional equation. 

For the GL(A) L— functions, the symmetric square has degree 10. Such high degree 
L— functions require a large number of coefficients to be computed accurately, so we were 
not able to obtain convincing results using the symmetric square. But for the specific case 
of Sp(4, Z) there is another option. We claim that we have found the spin L— functions as- 
sociated to Maass forms on Sp(4, Z). If that is true, then we can use the functional equation 
data and Dirichlet coefficients to compute the standard L— function, which has degree 5. 
(That the standard L— function can be reconstructed from the spin L— function can be seen 
from an examination of the explicit representations given in [13] .) We computed the degree 5 
L— functions in several cases and found that they satisfy the expected functional equation 
to high accuracy. We consider this to be a strong confirmation of the legitimacy of our 
calculations. Numerical data will be made available at |http: / / www. L MFDB . org/L / degree5 . 



7. Proof of Theorem 13.11 



In this section we describe the proof that the lighter shaded region in Figure 3.1 contains 



no degree 3 L— functions. The idea of the proof is to use the approximate functional equation 
method of Section [4] to form a system of linear equations in the Dirichlet coefficients, and 
then verify that the solution to that system violates the Ramanujan bound, when (Ai,A2) 
are in the shaded region. 

We will write out the details for the point (Ai, A2) = (8.4, 14.2). Note that this point is in 
our shaded region, but it is not in the region excluded by Miller's [20] argument using the 
explicit formula. 



In (4.10) we used the point (Ai, A 2 ) = (8.4, 14.2) to form one equation for the coefficients. 



Now we will make another equation. Choosing gs(s) = e lS//4 , and continuing to use so - „ 
we find 



g 3 (s ) 1 RHS(s 0) Q, «, A, e, {a n }, g 3 ) = (7.1) 
- 1078.21 - 2.38608 61 (2) - 256.449 6 2 (2) - 3.595 b x (3) + 17.50 3 6 2 (3) 
+ 1.0549 6i(4) - 0.1226 6 2 (4) + 4.84 • 10~ 6 61(8) + 0.000012 6 2 (8) + • • • . 

Recall that a(n) are the Dirichlet coefficients, and we write a(n) = b\(n) + i6 2 (n) with 



Since (4.8), (4.9), and (7.1) each equal Z{\ + i), setting them all equal gives two linear 
equations in the Dirichlet coefficients. We can solve for &i(2) and 6 2 (2) to obtain 

& x (2) = 0.1866 + 0.0137 6i(3) + 0.22 9 6 2 (3) + 0.0 1 2 2 6i(4) + 0.0380 6 2 (4) + ••• 

- 2.7 ■ 10" 7 61(8) - 2.5 ■ 10~ 6 62(8) + 2.3 ■ 10~ 13 6 X (14) - 1.6 ■ HT 12 6 2 (14) ■ ■ ■ . 

6 2 (2) = - 4.1854 - 0.0170 6i(3) + 0.0629 6 2 (3) + 0.0034 61 (4) - 0.0015 6 2 (4) + • • • 

- 2.4 ■ 10~ 8 61(8) + 1.4 • 10~ 7 62(8) + 1.3 ■ 10" 14 6i(14) + 6.1 ■ 10" 14 6 2 (14) 

(7-2) 
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Now we impose the Ramanujan bound (3.4) on all the coefficients on the right side of the 
above equation, giving 

h (2) = 0.1866 ±1.33 

6 2 (2) = -4.1854 ±0.301. (7.3) 

Note that the ± terms are rigorous bounds, not error estimates. Therefore 62(2) < —3.8842, 
provided all the other Dirichlet coefficients satisfy the Ramanujan bound. Since the Ra- 
manujan bound |a(2)| < 3 implies |&2(2)| < 3, we have a contradiction. The conclusion 
is that if (Ai,A2) are the parameters in the T-factors of an L— function, then not all the 
Dirichlet coefficients satisfy the Ramanujan bound. Thus, assuming the Ramanujan bound, 
there is no L— function with (Ai, A 2 ) = (8.4, 14.2). 

We eliminated (Ai, A 2 ) = (8.4, 14.2) by solving only two linear equations. Identical calcu- 



lations over a closely spaced grid of points give the shaded region in Figure 3.1 This proves 
Theorem 13.11 

It is possible to eliminate a larger region by solving a larger system, but we have not 
attempted to explore the limit of this method. 
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